arXiv: 1502.02673v2 [quant-ph] 18 Sep 2015 


Coherence and measnrement in quantum thermodynamics 

Philipp KammerlandeiQ 

Institute for Theoretical Physics, ETH Zurich, Wolfgang-Pauli-Strasse 27, 8093 Zurich, Switzerland. 

Janet Ander^ 

Department of Physics and Astronomy, University of Exeter, Stocker Road, EXf 4QL, United Kingdom. 

Thermodynamics is a highly successful macroscopic theory widely used across the natural sciences 
and for the construction of everyday devices, from car engines and fridges to power plants and solar 
cells. With thermodynamics predating quantum theory, research now aims to uncover the thermo¬ 
dynamic laws that govern finite size systems which may in addition host quantum effects. Here we 
identify information processing tasks, the so-called “projections”, that can only be formulated within 
the framework of quantum mechanics. We show that the physical realisation of such projections 
can come with a non-trivial thermodynamic work only for quantum states with coherences. This 
contrasts with information erasure, first investigated by Landauer, for which a thermodynamic work 
cost applies for classical and quantum erasure alike. Implications are far-reaching, adding a thermo¬ 
dynamic dimension to measurements performed in quantum thermodynamics experiments, and pro¬ 
viding key input for the construction of a future quantum thermodynamic framework. Repercussions 
are discussed for quantum work fluctuation relations and thermodynamic single-shot approaches. 


I. INTRODUCTION 

When Landauer argued in 1961 that any physical real¬ 
isation of erasure of information has a fundamental ther¬ 
modynamic work cost he irrevocably linked thermody¬ 
namics and information theory [T]. A practical conse¬ 
quence of this insight is that all computers must dissipate 
a minimal amount of heat in each irreversible computing 
step, a threshold that is becoming a concern with future 
computer chips entering atomic scales. The treatment of 
general quantum information processing tasks within the 
wider framework of quantum thermodynamics has only 
recently begun. Theoretical breakthroughs include the 
characterisation of the efficiency of quantum thermal en¬ 
gines nn and the extension of widely used classical non¬ 
equilibrium fluctuation theorems to the quantum regime 
011]. A new thermodynamic resource theory [7| has 
led to the discovery of a set of second laws that replaces 
the standard macroscopic second law for finite size sys¬ 
tems These results have substantially advanced 

our understanding of nanoscale thermodynamics, how¬ 
ever putting a finger on what is genuinely “quantum” 
in quantum thermodynamics has remained a challenge. 
Quantum mechanics differs from classical mechanics in 
at least three central aspects: the special nature of mea¬ 
surement, the possibility of a quantum system to be in 
a superposition and the existence of quantum correla¬ 
tions. The thermodynamic energy needed to perform a 
(selective) measurement has been investigated [H] and 
the total work for a closed thermodynamic measurement 
cycle explored m- The catalytic role of quantum su¬ 
perposition states when used in thermal operations has 
been uncovered HH and it has been shown that work can 
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FIG. 1: Thermodynamic setting. A system, depicted 
as a spin, interacts with a heat bath at temperature T, with 
which it exchanges heat, and with controlled energy sources, 
illustrated as coil and weight, with which it exchanges work. 
Work drawn from the system can be collected in a work stor¬ 
age system (weight) for future use. 


be drawn from quantum correlations [HE] in a ther¬ 
modynamic setting, see Fig. In particular, del Rio et 
al. m showed that contrary to Landauer’s principle, it 
is possible to extract work while performing erasure of a 
system’s state when the system is correlated to a mem¬ 
ory. This can occur if and only if the initial correlations 
imply a negative conditional entropy, a uniquely quan¬ 
tum feature. The thermodynamic process does however 
now require operation on degrees of freedom external to 
the system, i.e. the memory’s. 

Our motivation is here to shed light on the implica¬ 
tions of performing a measurement on a quantum state 
that has coherences. We will consider this task in the 
thermodynamic setting of Landauer’s erasure, involving 
a heat bath at fixed temperature T and operation on 
Af —>■ oo uncorrelated and identically prepared copies 
of the system (i.i.d. limit). This is of interest in the 
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context of the quantum Jarzynski equality, for example, 
and will also be central for experiments testing quantum 
thermodynamic predictions in the future. To tackle this 
question we define the information-theoretic “projection” 
p —t 77 ^ := for a given initial quantum state 

p and a complete set of mutually orthogonal projectors 
{n^}fc- Such state transformation can be seen as analo¬ 
gous to the state transfer of erasure, p —)■ | 0 ), to a blank 
state |0). Physically, this projection can be interpreted 
as the result of an unread, or unselective [15], measure¬ 
ment of an observable V that has eigenvector projectors 
{nHfc- In an unselective measurement the individual 
measurement outcomes are not recorded and only the 
statistics of outcomes is known. In the literature the 
implementation of unselective measurements is often not 
specified, although it is typically thought of as measur¬ 
ing individual outcomes, e.g. with a Stern-Gerlach ex¬ 
periment, see Fig. [ 2 ^, followed by mixing. The crux is 
that the information-theoretic projection p —>■ can 

be implemented in many physical ways. The associated 
thermodynamic heat and work will differ depending on 
how the projection was done and we will refer to the 
various realisations as “thermodynamic projection pro¬ 
cesses”. One possibility is decohering m the state in 
the so-called pointer basis, { 11 ^°“*®’^}^, a thermodynamic 
process where an environment removes coherences in an 
uncontrolled manner resulting in no associated work. In 
general it is possible to implement the state transfer in 
a finely controlled fashion achieving optimal thermody¬ 
namic heat and work values. 


II. MAIN RESULT 

Of particular importance in thermodynamics is the 
projection p —>■ rj^ of the system’s initial state p onto the 
set of energy eigenstates {^^}k of the system’s Hamil¬ 
tonian H — J2k^k ^k with Ek the energy eigenvalues. 
Here the state’s off-diagonals with respect to the energy 
eigenbasis are removed - a state transformation that is 
frequently employed in quantum thermodynamic deriva¬ 
tions and referred to as “dephasing” or “measuring the 
energy”. Our key observation is that there exists a ther¬ 
modynamic projection process realising this transforma¬ 
tion and allowing to draw from the quantum system a 
non-trivial optimal average work of 

{Wopt) = kBT {S{r,^)-S{p)). (1) 

Here T is the temperature of the heat bath with which 
the system is allowed to interact, see illustration Fig. 
ks is the Boltzmann constant and S is the von Neumann 
entropy. Crucially, this work is strictly positive for quan¬ 
tum states with coherences. Extending the key observa¬ 
tion to general projections p —>■ one finds that optimal 
thermodynamic projection processes can be implemented 
that allow to draw an average work of 

(VFopt) =kBT ( 5 ( 77 ^) - S{p)) - tr [H ( 77 ^ - p)] , (2) 


where an additional internal energy change term appears. 

The optimal work values stated in Eqs. Q and ([^ 
are valid for processes applied to classical and quantum 
states alike. While for a classical ensemble the entropy 
change, AS^ = S{rj^) — S{p), will be zero this is not 
so in the general quantum situation, where initial non¬ 
diagonal quantum states result in a strictly positive en¬ 
tropy change m We note that while the optimal work 
values are in principle attainable, practical implemen¬ 
tations may be suboptimal resulting in a reduced work 
gain or a higher work cost. The physical meaning of AS^ 
can be grasped by considering a lower bound [TB] on it, 
AS^ > I ||p — see AppendixHere d is the 

dimension of the system and || • II 2 denotes the Hilbert- 
Schmidt norm. The first factor quantifies the distance 
of the initial state from the fully mixed state, while the 
second factor, AA^, quantifies the angle between the di¬ 
agonal basis of p and the projection basis {^}k- These 
terms correspond to incoherent and coherent mixing con¬ 
tributions. The entropy change is non-trivially bounded 
only if the initial state is not an incoherent mixture with 
respect to that basis. The entropy bound is the largest 
for pure initial states whose basis is mutually unbiased 
with respect to {n^}fe. In this case the optimal entropy 
change is AS'’^ = kBT Ind. 

One may wonder where the work has gone to. There 
are two equivalent approaches to the accounting of work. 
In the present analysis the focus is on the work that the 
system exchanges, as done in statistical physics pm- 
I22j . In this approach it is often not explicitly mentioned 
where the work goes to, but the only place work can go 
to are the externally controlled energy sources. Similarly, 
the heat, i.e. the energy change minus the work, is es¬ 
tablished implicitly. For example, in the experimental 
realisation of classical Landauer erasure with a colloidal 
silica bead trapped in an optical tweezer [ 21 ], the dis¬ 
sipated heat of erasure was calculated by knowing the 
applied tilting forces and integrating over the bead’s dy¬ 
namics. The second approach is to collect work in a 
separate work storage system [23], as illustrated by the 
weight in Fig. [^ and detailed in Appendix [C] Both the 
implicit and the explicit treatment of work are equivalent 
in the sense that the results obtained in one approach can 
be translated into the other. 

The thermodynamic assumptions made to prove 
Eq. ([^ are congruent with current literature [DIMIS]; 
specifically they are: (TO) an isolated system is a system 
that only exchanges work and not heat; (Tl) the valid¬ 
ity of the first law relating the internal energy change, 
AU, of the system during a process to its average heat 
absorbed and work drawn, AU = — (IF); (T2) the 

validity of the second law relating the system’s entropy 
change to its average absorbed heat, kBT AS > (Q®''’’®), 
when interacting with a bath at temperature T, with 
equality attainable by an optimal process; (T3) the ther¬ 
modynamic entropy to be equal to the von Neumann en¬ 
tropy in equilibrium as well as out-of-equilibrium, ^th = 
^vN- In addition we make the following standard quan- 
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FIG. 2: Two physical realisations of a projection process, a, N identically prepared spin 1/2 particles in state 
pass a Stern-Gerlach magnet and a screen after which they emerge in either the spin-up or the spin-down beam. Recombining 
the two beams mixes the spins to the final state for N —>■ oo. Illustration of the spin example, showing the state evolution 
in b and the B-field evolution in c. The poles in the Blochsphere (b) are the energy eigenstates |eo) and |ei) that are aligned 
and anti-aligned with an externally applied B-field (indicated in blue in c), which initially is (black point in c). In the first 
step the Blochvector Sp (black arrow in b) of Emmy’s initial state p is rotated on the green-dashed circle to s/j (green arrow in 
b). The unitary rotation V reqnired for this step can be realised by applying a microwave pulse creating an additional B-field 
(indicated in orange in c) in the direction orthogonal to the plane of the green circle. At the end of the first step the pulse is 
turned off and the external B-field is adjusted to (green point in c). The second step shortens s/j to Sp (red arrow in b), 
the Blochvector of p (superscripts H have been omitted). The external B-field (blue in c) decreases slowly to (red point 
at t 2 in c). In the last step the B-field quickly returns to its initial value, (red point at ta in c), while the state remains 
p. The angle between the Blochvectors of p and p is indicated by 6. 


turn mechanics assumptions: (QO) an isolated system 
evolves unitarily; (QI) control of a quantum system in¬ 
cludes its coherences. Details of the proof are in Ap¬ 
pendix]^ We note that in the single-shot setting whole 
families of second laws apply O [9] that differ from (T2) 
stated above. However, in the limit of infinitely many in¬ 
dependent and identically prepared copies of the system 
these collapse to the standard second law, (T2), on the 
basis of which Eq. (§ is derived. 


From the information-theory point of view the pro¬ 
jections considered here constitute just one example of 
the larger class of trace-preserving completely positive 
(TPCP) maps characterising quantum dynamics. Of 
course, all TPCP maps can be interpreted thermody¬ 
namically with the assumptions stated above, resulting 
in an optimal average work given by a free energy differ¬ 
ence. Erasure is another such map whose study forged 
the link between information theory and thermodynam¬ 
ics. The benefit of discussing “projections” here lies in 
the insight that this focus provides: it uncovers that co¬ 
herences offer the potential to draw work making it a 
genuine and testable quantum thermodynamic feature. 
This work is non-trivial even when the thermodynamic 
process is operated on the system alone, not involving 
any side-information El stored in other degrees of free¬ 
dom. 


III. EXAMPLE 

To gain a detailed understanding of thermodynamic 
projection processes that give the optimal work stated in 
Eq. (|^ we now detail one such process for the example of 
a spin-1/2 particle (qubit), see illustration in Fig.]^ and 

as well as Appendix]^ This process consists of a uni¬ 
tary evolution, a quasi-static evolution and a quench [3S] , 
and it is optimal for any finite-dimensional quantum sys¬ 
tem as shown in Appendix]^ An experimentalist, Emmy, 
prepares the spin in a state p = a|0)(0| + (1 — a)|l)(l| 
(a > ^ w.l.o.g.) exposed to an external magnetic field 
B(o) -^bich she controls. The Hamiltonian associated 
with the system is H = —E|eo)(eo| + E|ei)(ei| where 
the energy difference between the aligned ground state, 
|eo), and anti-aligned excited state, |ei), is given by 
2E = \p\\B<-°'>\ with jl the spin‘s magnetic moment. 
Importantly, in general the spin state’s basis, {|0), |1)}, 
are superpositions with respect to the energy eigenba- 
sis, |0) = a* |eo) + /3* |ei) and |1) = j3 |eo) — a. |ei) with 
|ap -b |/3p = 1. For the optimal implementation of the 
projection p ^ p^ = X)fe=o,i |efc)(efc| P |efc)(efe| Emmy 
now proceeds with the following three steps. 

Firstly, she isolates the spin from the bath and mod¬ 
ifies external magnetic fields to induce a unitary rota¬ 
tion, V = |eo)(0| + |ei)(l|, of the spin into the energy 
basis. In nuclear magnetic resonance (NMR) [26l and 
pulsed electron spin resonance (ESR) experiments [57] 
such rotations are routinely realised by radio-frequency 
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and microwave pulses respectively, as evidenced by Rabi 
oscillations. The power, duration and phase of such a 
pulse would be chosen to generate the spin-rotation along 
the green circle until the desired unitary V is achieved. 
In the same step Emmy adjusts the strength of the ex¬ 
ternal B-field such that the spin state pi = V pV^ is 
Boltzmann-distributed at temperature T with respect to 
the energy gap of the Hamiltonian at the end of the step. 
In NMR or ESR the B-field magnitude is tuned 
quickly on the Ti timescale to achieve the desired energy 
gap. In the second step, Emmy wants to implement a 
quasi-static evolution of the spin that is now thermal. 
She brings the spin in contact with the heat bath at tem¬ 
perature T and quasi-statically adjusts the magnitude of 
the external B-field allowing the spin state to thermalise 
at all times. The final B-field, B^'^\ is chosen such that 
the final thermal state becomes . In ESR this step 
can be realised by changing the external B-field slowly 
on the Ti timescale so that the spin continuously equili¬ 
brates with its environment. Einally, Emmy isolates the 
spin from the environment and quickly changes the B- 
field to its original magnitude while the state remains 
. 

During Step 1 and 3 the system was isolated and 
the average work drawn is thus just the average energy 
change. During Step 2 the average work is the equi¬ 
librium free energy difference between the final and ini¬ 
tial thermal states at temperature T, see Appendix 
for details. In NMR/ESR the work contributions drawn 
from the spin system are done on the external B-field 
and the microwave mode. This could be detected by 
measuring the stimulated emission of photons in the mi¬ 
crowave mode or observing current changes induced by 
the spins dynamics [11I17]. The overall thermodynamic 
process has now brought the spin from a quantum state 
with coherences, p, into a state without coherences, 77 '^, 
while keeping the average energy of the spin constant. 
The net work drawn during the three steps adds up to 
(IT) = ksT (S(p^) — S{p)) showing the attainability of 
the optimum stated in Eq. Q for the spin-1/2 example. 
We note that Eq. 0 is also the maximal work that can 
be extracted from a qubit state p under any transfor¬ 
mation of the system that conserves its average energy, 
U := tr[H p], i.e. for qubits p^ is the optimal final state 
under this condition. 

We emphasise that this optimal implementation in¬ 
volves a finely tuned and controlled operation that relies 
on knowledge of the initial state p. This is akin to the 
situation considered in [14j where knowledge of the initial 
global state of system and memory is required for optimal 
erasure with side-information. It is important to distin¬ 
guish this situation from that of Maxwell demon’s who 
has access to knowledge of the individual micro-states 
that make up the ensemble state p^, and who uses 
it to beat the second law [28] . In the scenario considered 
here there is no knowledge of the individual micro-states 
and the process does not violate the second law, on 
the contrary, it is derived from it. 


IV. IMPLICATIONS 
A. Single-shot analysis 

The preceding discussion concerned the average work 
that can be drawn when operating on an ensemble of 
TV —)■ 00 independent spins. This scenario contrasts 
with the single shot situation considered in a number 
of recent publications [HUHIMIEDI- In particular, two 
major frameworks |29l I30j have recently been put for¬ 
ward to identify optimal single-shot work extraction and 
work cost of formation in the quantum setting. These 
frameworks rely on a resource theory approach |7] and 
make use of min- and max-relative entropies that orig¬ 
inate from one-shot information theory. The optimal 
work extraction schemes of these frameworks require non¬ 
diagonal states to be decohered first to become diagonal 
in the energy basis. This decoherence step is assumed 
to not have an associated single-shot work. However, 
the present analysis of energy basis projections showed 
that thermodynamic projection processes can yield pos¬ 
itive average work, see Eq. Q. Therefore one may ex¬ 
pect a positive work for removing coherences from a state 
p in the single-shot setting, too. Since our focus is the 
iV —>■ 00 limit we will not aim to construct the single-shot 
case. Nevertheless, to establish a notion of consistency 
between single-shot results [29l [30] and the average anal¬ 
ysis presented here we now separate the projection into 
a diagonal part that can be analysed in the single-shot 
framework and a non-diagonal part that can be anal¬ 
ysed in the average framework. One possible decompo¬ 
sition of p —>■ ? 7 ^ is the split in three steps each starting 

and ending with Hamiltonian H\ p pi \ A p^. 
Here pi is the rotated state defined above and = 
jg thermal state for the Hamil¬ 
tonian H at temperature T. We can now use a single-shot 
analysis [501 for Steps b and c that involve only states 
diagonal in the energy basis, giving a single-shot work 
contribution of /cbT In 2 (Dmin(pi 11 

see Appendix]^ Here and Dmax are the min- and 

max-relative quantum entropies, respectively. Taking the 
limit of iV —>■ 00 copies for Steps b and c and adding the 
average work contribution for the initial non-diagonal ro¬ 
tation a, (irfo)) = —tr[(pi — p)H], one indeed recovers 
the optimal average work as stated in Eq. 0 . After mak¬ 
ing public our results very recently a paper appeared [5T] 
that derives the work that can be extracted when remov¬ 
ing coherences in a single-shot setting. These results are 
in agreement with Eq. 0 and reinforce the above con¬ 
clusion that coherences are a fundamental feature distin¬ 
guishing quantum from classical thermodynamics. 


B. Quantum fluctuation relations 

The key observation was that thermodynamic projec¬ 
tion processes can have a non-trivial work and heat. An- 
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other instance where this has interesting repercussions is 
the quantum Jarzynski equality mi. This is a gener¬ 
alisation of the prominent classical fluctuation relation 
valid for general non-equilibrium processes, which has 
been used to measure the equilibrium free energy sur¬ 
face inside bio-molecules by performing non-equilibrium 
pulling experiments |19j . The quantum version has re¬ 
cently been tested for the first time in a nuclear mag¬ 
netic resonance experiment |2b) . The quantum Jarzynski 
relation, links the fluctuat¬ 

ing work, W, drawn from a system in individual runs of 
the same non-equilibrium process, with the free energy 
difference, AF, of the thermal states of the final and 
initial Hamiltonian, see Appendix m In its derivation 
a system initially in a thermal state pq with respect to 
Hamiltonian at temperature T is first measured in 
the energy basis of The Hamiltonian is then varied 

in time ending in generating a unitary evolution, 

V, of the system, see Fig. [3}r. A second measurement, 
in the energy basis of is then performed to estab¬ 

lish the final fluctuating energy. For each run the differ¬ 
ence of the two measured energies has been associated 
with the fluctuating work |5], AE = —W. The experi¬ 
ment is repeated, each time producing a fluctuating work 
value. On average the work extracted from the system 
during the quantum non-equilibrium process turns out 
to be = U{po) — U{pr) where pr '■= VpoV^ 

is the ensemble’s state after the unitary evolution, and 
similarly the average exponentiated work is calculated. 
The above identification W := —AE was made assum¬ 
ing that the system undergoes a unitary process with no 
heat dissipation. However, the need to acquire knowl¬ 
edge of the system’s final energies requires the second 
measurement. The ensemble state is thus further altered 
from Pt to Pr, the state pr with any coherences in the 
energy basis of removed. This step is not unitary 
- during the projection pr —>■ 77 ,-the system may absorb 
heat, (Q®'*’®), indicated in Fig. |^, whose value depends 
on how the process is conducted. Thus, while the energy 
difference for the projection is zero, U{pr) — U{pr) = 0 , 
for states Pr with coherences the entropy difference is 
not trivial, S{pr) - S{pr) = {QopD/iksT) > 0. This 
implies that in an experimental implementation of the 
Jarzynski relation the work done by the system on av¬ 
erage can be more than previously thought, (Wopt) = 
+kBT {S{pr) — S{pr))- We conclude that the 
suitability of identifying W = —AE, and hence the valid¬ 
ity of the quantum Jarzynski work relation, depends on 
the details of the physical process that implements the 
second measurement. This conclusion is not at odds with 
previous experiments | 2 tij which showed nature’s agree¬ 
ment with = g-AF/(fci 5 T)^ involving the 

average of the exponentiated measured fluctuating en¬ 
ergy. 
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FIG. 3: Dynamical steps in a quantum fluctuation ex¬ 
periment. a, The quantum Jarzynski relation is described 
as characterising the non-equilibrium work of processes that 
start in a thermal state po and evolve unitarily (V), driven by 
a changing Hamiltonian, reaching the final state pr at time r. 
This unitary process has no heat contribution, b, Illustration 
of three steps that are assumed in mathematical derivations of 
the quantum Jarzynski relation mi]: initial energy measure¬ 
ment of indicated by Mq, unitary evolution, and final 
energy measurement of JjA) indicated by Mr- The ensemble 
state evolves here from po to pr and then to pr, the state pr 
with its coherences removed. The observed average energy 
difference AU = U{pr) — U{po) encompasses both, the uni¬ 
tary process and the second projection process, and can in 
general contain a heat contribution in contrast to a. 


C. Correlated systems 


It is insightful to extend the thermodynamic analysis 
of projections to correlated systems. An experimenter 
may have access not only to the system S but also the 
auxiliary systems A with which S is correlated [14) . She 
can then perform a global operation, p^"^ —)■ that 

implements a projection {H^j^ locally on the system S, 
i.e. p^ := trA[p®^] ^ p^ := tr a[p^'^] = '^'k , 

while leaving the reduced state of the auxiliary system 
unchanged, i.e. p^ := trs[p‘®^] = tr 5 [ 77 ‘^^]. By doing 
so the experimenter can optimally draw the overall work 
(Wopt) = ksT AS^ — At/^, where AS^ is the entropy 
change for the state of system-bauxiliary and AC/^ is still 
the energy change of the system alone. This quantity can 
be re-written as the sum of two terms: (Wopt), the ex¬ 
tractable work when operating on the system alone given 
in Eq. ([^, and {A : S), a positive term quantifying the 
quantum correlations between S and A, see Appendix [H| 
The latter contribution was previously identified in an in¬ 
spiring paper by Zurek m- It depends on the choice of 
projectors and is related to, but broader than, quantum 
discord [32] which is optimised over all possible projec¬ 
tors. This means that even states of system and auxiliary 
that can be considered classically correlated (i.e. no dis¬ 
cord) provide an advantage for drawing work contrasting 
with the erasure process where this only occurs for highly 
entangled states [14]. The gap between these two sets of 
correlated states is an intriguing fact and calls for further 
exploration of the link between thermodynamics and in¬ 
formation theory in the quantum regime. 
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V. CONCLUSIONS 

To conclude, erasure is not the only irreversible infor¬ 
mation processing task - in the quantum regime a second 
fundamental process exists that mirrors Landauer’s era¬ 
sure. In contrast to the minimum heat limit of erasure, 
thermodynamic projection processes have a maximum 
work limit. While the former is non-zero for the erasure 
of classical and quantum bits, optimal thermodynamic 
projection processes have a non-zero work only when ap¬ 
plied to quantum states with coherences. The optimal 
average work stated in Eqs. 0 and 0 constitutes an 
experimentally accessible quantum thermodynamic pre¬ 
diction. Future experiments testing this optimal work 
may be pursued with current setups, for instance with 
NMR/ESR techniques [551 EZ] or single atoms [551 IMj . 
and promise to be accessible with other platforms en¬ 
tering the quantum regime, such as single electron boxes 
[22] . Experiments will be limited by practical constraints, 
such as achieving a quasistatic process and obtaining the 
maximum work for pure states which may require, for 
instance, very large B-fields. 

The derivation of the optimal work value is mathe¬ 
matically straightforward, just like that of Landauer’s 
principle. The result’s significance is that it opens new 
avenues of thought and provides key input for the con¬ 
struction of a future quantum thermodynamic frame¬ 
work. For example, the developed approach opens the 
door to investigate the connection between microscopic 
statistical physics and macroscopic thermodynamics in 
the quantum regime. While it is straightforward to iden¬ 
tify the thermodynamic work of quantum processes in¬ 
volving macroscopic ensembles, what is needed is a mi¬ 
croscopic concept of work that when averaged, gives the 
correct macroscopic work. The microscopic work concept 
should be valid for general (open) quantum processes and 
quantum states (including coherences), and only require 
access to properties of the system. While single-shot ap¬ 
proaches have discarded coherences [551130] , fluctuating 
work approaches cannot be applied directly to a system 
undergoing open quantum evolution [50] . 

The observation is also important from the experimen¬ 
tal perspective as testing quantum thermodynamic pre¬ 
dictions will involve measurement - a projection process. 
We have argued that measurements, such as those re¬ 
quired in establishing the Jarzynski equality, are not nec¬ 
essarily thermodynamically neutral. Indeed, they can be 
implemented in different physical ways and in general 
play an active role in thermodynamics, contributing a 
non-zero average heat and work. This new perspective 
gives physical meaning to the change of entropy in the 
debated quantum measurement process - it provides a 
capacity to draw work. Specihcally, work can be drawn 
when coherences of a state are removed during an unse- 
lective measurement. 

Finally, it is apparent that optimal thermodynamic 
projection processes require use of knowledge of the ini¬ 
tial state p, i.e. its basis and eigenvalues. One may 


be inclined to exclude use of such knowledge, partic¬ 
ularly when considering projections in the context of 
measurement which is often associated with the acqui¬ 
sition of knowledge. Such restriction would necessar¬ 
ily affect the set of assumptions (T0-T3, QO-QI) in the 
quantum regime. These could be changed, for exam¬ 
ple, to the second law not being possible to saturate (cf. 
T2) or to choosing a new quantum non-equilibrium en¬ 
tropy that only considers the state’s diagonal entries (cf. 
T3). The latter would mean a departure from standard 
quantum information theory where entropies are basis- 
independent. Thus whichever approach one takes - not 
making or making a restriction - quantum coherences will 
contribute a new dimension to thermodynamics. They 
either lead to non-classical work extraction or they alter 
the link between information theory and thermodynam¬ 
ics in the quantum regime. The line drawn here between 
the assumptions (T0-T3, QO-Ql) and results (Eqs. 0 
and ([^) establishes a frame for this possibility to be in¬ 
vestigated. 
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Appendix A: Proof of Eq. ([^ 

Using the first law (TI) the average work drawn in 
a thermodynamic projection process p —>■ rj^ is simply 
{W) = (Q®'*’®) — AC/^, where A[/^ is the average energy 
change for that process. Relating the average heat ab¬ 
sorbed by the system during the process to its entropy 
change one then obtains {W) < ksT AS^ (T2). Here 
AS'P is the difference of von Neumann entropies of the 
system’s state before and after the projection (T3). The 
average work drawn is thus {W) < ksT AS^ — AC/^, 
where the entropy change is non-negative and the energy 
change can be either positive or negative. The stated 
optimal work, (Wopt)? is achieved when the inequality is 
saturated by an optimal process (T2) the implementa¬ 
tion of which may require knowledge of the initial state 
and control of coherences (Ql). In the special case of a 
projection onto the energy eigenbasis the inter¬ 

nal energy change is zero, AU^ = 0, and one obtains 
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Eq. 0. 

Appendix B: Spin example 

We here detail the three steps of the optimal thermo¬ 
dynamic projection process of the spin system discussed 
in the main text. Emmy starts with a spin in state 

p = a|0)(0| + (l-a)|l)(l|, (Bl) 

with Blochvector 

Sp := trjyoa] = (2a — 1) 6, (B2) 


Since geometrically the mapping p —>■ r]^ is a projec¬ 
tion of Sp onto the vertical axis in the Blochsphere, the 
length of the final Blochvector, s^, is shorter than the ini¬ 
tial Blochvector, Sp. This shortening is associated with 
an entropy increase When describing the process in 
the following we assume that p > ^ in accordance with 
the illustration in Fig. 0. At the end of this section we 
come back to the case p < 

Emmy proceeds with three steps made up of quantum 
thermodynamic primitives with known work and heat 
contributions [IS], Fig. 0: 

{p,H) ^ ^ A 


where a is the vector of the three Pauli matrices, cri,(T 2 
and (73 and 0 = tr[|0)(0|(T] is the unit vector in the 
Blochsphere pointing from the origin to the state |0), 
see Fig. We assume without loss of generality that 
a > i. If this was not the case, the labels |0)(0| and |1)(1| 
should be interchanged. The spin’s initial Hamiltonian is 
given hy H = —E (H^ —Hf^), where = \ek){ek\ with 
fc = 0,1 are the rank-1 projectors onto the two energy 
eigenstates and E > 0. This Hamiltonian arises when 
the spin is exposed to an external magnetic field 
The energy separation of the aligned ground state, |eo), 
and anti-aligned excited state, |ei), is 2E = 2 \fl\ 
where p is the magnetic moment of the spin. A general 
initial state p is not diagonal in the basis {|eo), |ei)}, in 
other words the spin’s eigenstates are superpositions with 
respect to the energy eigenbasis, |0) = a* |eo) -|- /3* |ei) 
and |1) = ,5 |eo) — a |ei) with \a\^ + |/3p = 1. The spin’s 
Blochvector, 6, is then not parallel to the B-field, 
Emmy wants to obtain the state where the coherences 
with respect to the energy basis {|eo),|ei)} have been 
removed, 


nfpBf 

fc=0,l 

= ^ tr[nf p] Hf = p Ho^ + (1 - p) Hf, 

fc=0,l 


(B3) 


where p = tr[n^ p] is the probability for obtaining out¬ 
come I eg) in a measurement of H. Here r]^ has the same 
average energy as the initial state p, 


tr[p^H]= ^ trjnfpHfiJ] 

k^0,l 

= i/nf ] = tr[p7J], 

k^0,l 


(B4) 


where we used orthonormality and completeness of the 
projectors {n^,n(^}. The Blochvector of the final state 
is defined as 


s,, := tr[p^ a] = (2p - 1) eg (B5) 

where eg = tr[|eg)(eg| ct] is the unit vector in the 
Blochsphere pointing from the origin to the state jeg). 


In the first step, (p, H) —!-)■ (pi, Emmy isolates the 

spin from the bath and rotates the B-field such that the 
variation of the field induces a unitary transformation of 
the spin into the energy eigenbasis, with unitary V = 
|eg)(0| -I- |ei)(l|. The state after this step is 

pi = VpV^ = aUl^+ {l-a)U^. (B6) 

The B-field after this step, B^^\ is chosen such that the 
new Hamiltonian E[^^'> = —E^'^'> (H^ — Hf^) has eigen¬ 
values = IpI ■'^here ks is the 

Boltzmann constant and T is the temperature of the heat 
bath that Emmy will use in the next step. This choice 
of the B-field makes the state pi a thermal state with 

respect to at temperature T, i.e. pi = — with 

= tr[e“^'^' '] and inverse temperature jS = 

Since the system was isolated in the first step no heat 
exchange was possible and the entire average energy 
change of the system is drawn from the system as work 
(IH(i)) = -tr[pii7(i) — pH\. Physical constraints may 
make this process difficult to realise, for instance, pure 
initial states would require a B-field, of infinite mag¬ 
nitude because thermal states at any finite temperature 
are only pure if the energy gap is infinite. In this case 
there is a trade-off between the maximal magnitude the 
B-field can reach and the precision with which the pro¬ 
cess is carried out. In the following we assume that the 
maximal B-field is large enough to make the error in the 
precision negligibly small. 

In the second step, ^ {r]^,E[^^^), Emmy 

brings the spin in contact with the bath at temperature 
T, not affecting the spin’s state as it is already ther¬ 
mal. She then quasi-statically decreases the magnitude 
of the B-field, while keeping the system in contact with 
the bath at all times, such that the final Hamiltonian 
is E[^‘^'> = —E^‘^'> (n^ — nf^) where the B-field is chosen 
such that E^‘^'> = In where p is the probability 
of measuring —E in the initial state, p. The quasi-static 
evolution means that the system is thermalised at all 
times, arriving in the final state 

p-/3r/'"> 

—^=pno^ + (i-p)nf = 77" (B 7 ) 
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which is thermal with respect to where 

= tr[e“^^''^*]. This state is exactly rj^, 

the desired final state after the projection. The 
quasi-static process considered here has a known 
average work given by the free energy difference 

^ m - fW(pi)) 

where F^\rj^) = tT[r]^ — ksTS^rj^) and 

F^^\pi)=tT[p,H(^'>]-kBTSipi) are standard thermal 
equilibrium free energies and S is the von Neumann 
entropy defined by S{p^) = —tT[p^ log 77 ^] and likewise 
for pi. 

Finally, in the third step, {p^,F[^'^'>) —^ ( 77 ^,i7), 
Emmy isolates the spin from the bath and changes the 
energy levels of the Hamiltonian such that it becomes the 
initial Hamiltonian F[ again. This step is done quickly 
so that the state of the spin does not change. Because 
the system is isolated the energy change in this step is 
entirely due to work = —tr[ 77 ^ {H — F[^‘^'>)]. In 

total, this thermodynamic process has brought the spin 
from the quantum state (p, H) to the state ( 77 ^, H) while 
not changing the energy of the spin, tr[(p — F[] = 0. 

The overall average work drawn from the spin is 

(IF) = (IF^i)) + (IF(2)) -h (IF(3)) 

= —tr[pi + tr[pH] — + ksT S{p^) 

+ tr[pi - keT S{pi) - tT[p^ {H - 

= kBT{S{p^)-S{p^)) 

= kBTiSip^)-Sip)) = {Wopt). 

showing the optimality of the three step process for the 
spin example, cf. Eq. Q. 

The above example assumed p > Suppose now that 
the probability to find the final state 77 ^ in the ground 
state I Co) with respect to the Hamiltonian was smaller 
than to find it in the excited state |ei), i.e. p < ^. 
Proceeding through the three steps described one finds 
that the mathematics is exactly the same. In particular, 
after Step 2 77 ^ is a thermal state with respect to Ff^^^ 
at inverse temperature /3. The only difference occurs 
in the interpretation as for the Hamiltonian the 

ground state is |ei) because In < 0 is 

negative. This is feasible by making the B-field 
negative, thus swapping the ground and the excited 
state. Consequently the analysis above and the resulting 
expression of the total extracted work remain the same. 

The work extracted in the individual steps of the 
thermodynamic projection process can be either positive 
or negative, depending on the initial state p, the Hamil¬ 
tonian iJ and the temperature T of the heat bath. Their 
sum, (IF), is strictly positive whenever the initial state 
was not diagonal in the energy eigenbasis, a consequence 
of the entropy increase m from p to 77 ^. On the other 
hand for classical states - all diagonal in the energy basis 
- the optimal work for such a projection is always zero. 


Appendix extends the optimality proof of the above 
three step process to the general finite-dimensional case. 

A note on optimal work extraction at constant 
average energy. Assume we are given an initial state 
p and a non-degenerate Hamiltonian FI for a quantum 
system. The goal is to find the maximal work that can 
be obtained in a thermodynamic process that involves a 
heat bath at temperature T under the restriction that 
the average energy of the system after the process is the 
same as it was before the process, U := tr[piJ]. Using 
Eq. ([^ together with the condition that internal energy 
does not change this amounts to finding the maximum 
over the set of states tr with tr[CT H] = U, 

(l^max) _ niaxfcBr(^(a) - S^p)) 

(B8) 

= kB T{ma.xS{a) — S{p)). 

<7 

It is well-known that at a fixed expectation value of an 
observable H the Gibbs states a\ = e“'^^/tr[e“''^] are 
the states of maximal entropy gniEi]. Here the parame¬ 
ter A has to be chosen such that the energy of the Gibbs 
state matches U - therefore there is only one ax*, with A* 
such that tr[cT;,* H] = U, that gives the maximum here. 
The maximum entropy is then 

<S'(crA*) =-trlcTA* IntTA*] =-tr[(TA*(-A*H - lntr[e“^ ^])] 

= X*U + \ntr[e-^‘^] 

(B9) 

and the maximum average work that can be extracted 
from p at fixed average energy U is then 

(IF“ax) = kB T{X* U + Intrje-^'-f^] - S{p)), (BIO) 

For the special case that the system is a qubit (two- 
dimensional) the optimum Gibbs state for work ex¬ 
traction ax* is identical to the projected state = 
and the maximal work that can be 
drawn from a system starting in state p, while keep¬ 
ing its average energy fixed, is (IFjpt) in Eq. 0. To 
see this we expand H = H^^i + H^^i and U = 
p£;(o)_|_(i_p) 7 ^(i) -vvithp := tr[pn(°i]. Now here A* must 
be chosen such that ax* = e~^ '^/tr[e“''* = pH^^i + 

(I—p)n^^i, i.e. p = +6“"^*-®'^'), so 

that ax* has just the right energy tr[(TA* H] = U. On 
the other hand the projection state has the same expan¬ 
sion, p^ = pH^^i + (1 ~ p) H^^i = ax*- We note that 
this coincidence is not true for higher dimensional sys¬ 
tems where the energy-projected state p^ will in general 
have a non-monotonous, non-canonical distribution in its 
energy eigenbasis, while ax* must be Gibbs-distributed. 

Gonsidering the illustration in Fig. the qubit states 
a fulfilling the condition tr[cr H] are located on the plane 
which contains p and is perpendicular to the |eo)-|ei)- 
axis. On the other hand, in the Bloch picture a state 
has higher entropy the closer it is to the center of the 
sphere. Hence, the optimal hnal state when extracting 
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work from p while conserving the average energy of the SBW have to conserve the energy exactly, not only on 

system is the state p projected to the |eo)-|ei)-axis, i.e. average, which amounts to = 0. 


Appendix C: Work storage system 

In the previous section it was stated that work can 
be drawn from a quantum system when undergoing a 
thermodynamic projection process. But where has the 
work gone to? 

There are two approaches of accounting for work that 
are mirror images to each other. One approach 
[Ml EH EH Sa |43] focusses on the work that the system 
exchanges, as described above. Here it is often not ex¬ 
plicitly mentioned where the work goes to, but the only 
place it can go to are the externally controlled energy 
sources, see Fig. Another way of accounting is to ex¬ 
plicitly introduce a work system to store the work drawn 
[Ml [50] . One way of doing so in an average scenario is 
to introduce [M] a ‘suspended weight on a string’, de¬ 
scribed by a quantum system W, that could be raised or 
lowered to store work or draw work from it. Specifically, 
the Hamiltonian of the work storage system is defined 
as JI^ = mg X, representing the energy of a weight of 
mass m in the gravitational field with acceleration g at 
height X. In addition, an explicit thermal bath B is intro¬ 
duced gniiiT] consisting of a separate quantum system 
in a thermal (or Gibbs) state r^. Both, the explicit work 
storage system and the heat bath are illustrated in Fig.j^ 
In the latter approach the total system starts in a prod¬ 
uct state of system S (e.g. spin), bath H, and weight W, 
pSBW _ pS ^ ^ ^ which together undergo average 

energy conserving unitary evolution with V: 

pSBW ^ ^SBW pSBWyt^ (Cl) 

The assumption is that the total Hamiltonian is the sum 
of local terms, -|- . The average 

energy conservation constraint then reads tr[ — 

pSBW-^ jjSBWj _ Q |.]^g average work extracted to 
the work storage system, (VF), is identified with 

-p^^'^)H'^] = {W). (C2) 

Both the implicit and the explicit treatment of work are 
equivalent in the sense that the results obtained in one 
language can be translated in the other and vice versa. 
In particular, the implicit description used in this text 
[25| has an equivalent explicit formulation [23] . 

In the next section we will discuss single-shot ex¬ 
tractable work in a projection process. One possibility to 
define work in this context is to chose the explicit work 
storage system as a ‘work qubit’ with a specific energy 
gap which has to be in a pure energy eigenstate before 
and after the protocol [M] • This way it is guaranteed that 
full knowledge about its state is present at all times and 
the work is stored in an ordered form. In this scenario 
the allowed unitary operations V on the whole system 


Appendix D: General three step process 

It is straightforward to generalise the proof of optimal¬ 
ity from the two-dimensional spin-I/2 example to ther¬ 
modynamic projection processes in dimension d. Again 
the projectors {H^jfc map onto the energy eigenspaces 
of the Hamiltonian, H = where 

k = 1,..., d, are the energy eigenvalues. A general ini¬ 
tial state can be written as p = \ j){j\ where 

Oj > 0 are probabilities, rank- 

1 projectors on the corresponding eigenvectors |j), and 
j = I,..., d. A unitary operation, V, is now chosen such 
that it brings the initial configuration (p, H) into the new 
diagonal and thermal configuration (pi = V pV\H^^'>) 
where pi = Ofe and ■ The new 

energy eigenvalues, are adjusted such that the prob¬ 
abilities Ofc are thermally distributed with respect to 
for the bath temperature T. Adjusting the Hamiltonian 
eigenvalues while letting the state thermalise at all times 
now results in a isothermal quasi-static operation from 
(pi,7d(i)) to Hf). 

Here the new energy eigenvalues, Ej, , are chosen to 
be thermal (at T) for the state’s probabilities which are 
given by = tr[pn(^]. Finally, a quench brings the 
thermal configuration (rj^, quickly into the non¬ 

equilibrium state {g^,H). The average work for this 
overall process is (IT) = where (kF^^^) = 

—tr[id^^) pi — Hp] and (IF*-^^) = —tr[id rj^ — rj^] be¬ 

cause the first and third steps are unitary (QO-I-TO). The 
quasistatic step’s work is (MHHj (kF^^^) = —f!^^ + 
where F^'^ = tr[id(^^ Pi] —fcs T S'(pi) is the thermal equi¬ 
librium free energy for Hamiltonian and similarly, 

f!^'^ = tr[ij(^) rj^] — kBT S{ri^). Summing up and using 
tr[id (p — T]^)] = 0, one obtains 

(IF) = kB nsip^) - S{p)) = (IFopt), (Dl) 
concluding the optimality proof of the process sequence. 


Appendix E: Lower bound on entropy change 

The entropy change during a projection with projec¬ 
tors {H^ = \<pk){(l>k\}k can be lower bounded. In the fol¬ 
lowing, IIHII 2 = y^tilB^B] denotes the Hilbert-Schmidt 
norm of a linear operator B acting on a d-dimensional 
Hilbert space describing the quantum system of interest. 
The lower bound reads |3S] 


AS-^ := 5(p^) - S{p) > ^ 




(El) 
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Here, S is the von Neumann entropy, p the initial state 
and if’ = J2k the final state after the projection 

process. Furthermore, AA^ is the second smallest eigen¬ 
value of the matrix 1 — where M is the doubly 

stochastic matrix given by the entries Mki = \{(t>k\l)\‘^ 
and {| i® the eigenbasis of the initial state p. 

Cons ider ing the two main terms on the right hand side 
of Eq. ( |E1| separately, ||p — l/d||2 and AA^, it becomes 
apparent that they characterise different properties of the 
initial state. The first term measures the distance of p to 
the fully mixed state, 1/d, and quantifies the purity of 
p. It is maximal for all pure initial states and zero if and 
only if p = 1/d. In the special case of a spin-1/2 system 
it can be directly related to the length of the Bloch vector 
describing p in the Bloch representation, a link that will 
be established below. The second term, is related 

to the overlap of the eigenbasis of p, {| Oli) the pro¬ 
jective basis, {\4>k)}k- It is zero if they are the same and 
maximal if they are mutually unbiased [31 [37]. This can 
be seen as follows: if the two bases are the same, then 
the matrix M is a permutation and consequently M'^M 
is the identity. In this case, 1 — M’^ M is the zero matrix 
and thus AA^ = 0. If {\'pk)}k and {| l)}i are mutually 
unbiased, i.e. if they fulfil |(</fc|0P = 1/*^ fo'' all k,l, 
the matrix M and thus also M'^M is a rank-1 projector 
onto the space spanned by the vector (1,..., 1)^. Hence, 
1 — M'^M has eigenvalues {0,One finds that 
the second largest eigenvalue is = 1, which is also 
the maximal eigenvalue the matrix 1 — M'^M can have 

m- 

In the special case of the spin-1/2 system shown in 
Fig. a , the bound reads AS^ [p) > | |sp|^ sin^ d, where 
Sp is me Bloch vector of the initial state and 6 is the angle 
between the eigenbasis of p, {|0), |1)}, and the projective 
energy basis, {|eo), |ei)}. Let p = a |0)(0| -b (1 - a) |1)(1| 
be the initial state of the qubit. Furthermore, let rj^ = 
p|eo)(eo| + (1 — p)|ei)(ei| be the hnal state after the 
energy projection, where p = tr[ |eo)(eo| p] is the proba¬ 
bility to obtain |eo). As argued in Appendix [b] w.l.o.g. 
we can assume that a > |, p > | ■ In the Bloch represen¬ 
tation one can write p = i(l-l-s-CT) and = ^{1+t-a). 
Here we used a different notation for the Bloch vectors 
of p, s := Sp, and t := Sp, for readability. The Pauli 
matrices are self-adjoint and fulfil tr[aiaj] = 2dij. Hence 
we find 


Therefore 

|(eo|0)|^=tr[|0)(0||eo)(eo|] 


1 

_ 1 

“ 4 


, Si t., 

1 + TVjCTi + —^(Ti 


Sitj 


l®1 


2-b0-b0-b2 


\i\ * |s1|tl 

s ■ t 


(Jid j 


(E3) 


= + 


where summations over double indeces are assumed and 
6 is the angle between the two Bloch vectors of p and 
p^, see Fig. ir. By using orthonormality of both bases 
{|Z)}; and {\ek)}k it is shown that in the qubit case the 
matrices M and M'^M have the form 


= /^l(eo|0)P |(eo|l)P\ 1 A-bcosd l-cosd\ 

Vi(eii0)i^ i(eiil)py 2Vl-cosd l-bcosdj 

(E4) 


and 


M^M 


1 /l -b cos^ 9 1 — cos^ 9\ 

2 \1 — cos^ d 1 -b cos^ d J ■ 


(E5) 


Computing the eigenvalues {Ai, A2} of M'^M yields Ai = 
1 and A2 = cos^ d. Therefore AA^ = 1—A2 = 1—cos^ d = 
sin^ d. In total, this concludes the proof of the bound for 
the special case of a spin-1/2 system when projected in 
the energy eigenbasis {nf}fe = {|eo),|ei)}, 

AS^{p) > i |sp|^ sin^d. (E6) 

To further illustrate the bound consider the special 
case when the initial state p is pure and its eigenbasis 
mutually unbiased with respect to the energy eigenbasis, 
{|eo), |ei)}. In this case the final state after the projec¬ 
tion, p^, is maximally mixed and we hnd 

AS^{p) = S{r]^) - S{p) = ln2 - 0 « 0.69 . (E7) 

Here, the lower bound is equal to ^ =0.25 because |AI = 
1 for a pure state p and sin^ d = 1 for mutually unbiased 
bases. Thus in this example the bound is not particularly 
tight. 


Appendix F: Single-shot analysis 


P 


1 

2 


'1- 

(l^ A] 

— 

= tr 


7:5 • cr 

-S • cr I 

2 

2 


J 

V2 J\ 


^ ^ 1 1 
= 4'S*Sitr[cr,crj] = 

bi=i 


(E2) 


where | • | is the Euclidean metric in This proves the 
form of the first factor in the bound. Eor the factor AA^ 
notice that by assumption a > 5, p > \ and thus we can 

write |eo)(eo| = 5(1 -b ||y ct) and |0)(0| = |(1 + ||| • A- 


Instead of performing a thermodynamic process on an 
ensemble of N identical and independent copies one can 
consider a single run of the process. Two major recent 
frameworks [311130] have been developed to describe the 
optimal work that can be drawn from a system in a sin¬ 
gle run. The proposal by Aberg [H], involves changes 
of the Hamiltonian and identifies work with the deter¬ 
ministic energy change of the system when undergoing a 
unitary process. The proposal by Horodecki-Oppenheim 
[30) . is formulated in terms of thermal operations [44) . 
where work is associated with raising a two-level system. 
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a decohering 


iP:H) 

b basis rotation 

thermai operations 


{P,H) 

II 

iPuH) 



c 

quantum thermai operations 


{P,H) 



II 


state pi- As discussed, it cannot be treated with the 
single-shot framework [29l [30] but it is possible to asso¬ 
ciate an average extracted work with this unitary process, 
= tr[(p —pi) H]. A single-shot analysis according 
to Horodecki-Oppenheim ISO] can then be performed for 
the diagonal steps b and c. This is possible because the 
steps go via the thermal state . Step b brings pi to 
and allows the extraction of the single-shot work [50] 

M/W =fcBrin2i?A„(pi||r^), (FI) 


FIG. 4: Illustration of possible quantum thermodynamic 
processes that transform (p^H) to a, Decohering 

the state in the energy basis extracts no work, b, To perform a 
consistency check between the average and single-shot results it is 
possible to split the process into a basis rotation to (pi,H) with 
unknown single-shot work, but known average work, and two ther¬ 
mal operations that pass through the thermal state and 

are treatable in the single-shot framework Eq]. c, General quan¬ 
tum thermodynamic processes could allow coherences and need not 
pass through intermediate fixed states. 


called the ‘work qubit’, with energy gap W deterministi- 
cally from the ground to the excited state. 


However, when attempting to apply these two frame¬ 
works to find the single-shot work for the energy pro¬ 
jections p rj^ captured by Eq. Q one encounters an 
obstacle: both frameworks only apply to processes be¬ 
tween initial and final states that are classical, i.e. states 
that are diagonal in the energy basis. Aberg discusses co¬ 
herences in a separate framework |12j . which does how¬ 
ever not cover single-shot work extraction and only fo¬ 
cusses on average quantities, similar to those in other 
references US] US]. Horodecki-Oppenheim suggest that 
quantum states with coherences with respect to the en¬ 
ergy eigenbasis are first decohered before applying the 
single-shot protocol. As discussed, apart from decoher¬ 
ing there are other thermodynamic projection processes 
that map the initial state with coherences, p, to the fi¬ 
nal state P^k without coherences, where 

are the projectors on the energy eigenstates of the 
Hamiltonian, H. Eq. shows that the average work 
extracted in an optimal thermodynamic projection pro¬ 
cess is strictly positive while the decoherence process has 
zero work. Therefore one may expect a positive optimal 
work for projections also in the single-shot setting, with 
decohering a suboptimal choice, see Fig. [^ 


Since our focus here is the V —> oo limit we will not 
aim to construct the single-shot case. Instead, to estab¬ 
lish a notion of consistency between the average analysis 
and previous single-shot work results we consider the se¬ 
quence {p,H) in 

which the Hamiltonian before and after each step are the 
same and pi is the rotated state defined above, Eq. (B61. 
Here ^ ^ is the thermal state for Hamiltonian 

H at inverse temperature j3 = and Z = tr[e“^'^] 
is the partition function. Step a of this sequence ro¬ 
tates the initial non-diagonal state p to the diagonal 


where is the smooth min-relative entropy [45] and 
e > 0 is the allowed failure probability of the process. 
Similarly, in Step c the final state is formed from the 
thermal state by applying a protocol that costs work. 
This work is [5D] 

IT(=) = -fcBTln2i?^,,(ry^||r^), (F2) 

where is the smooth max-relative entropy |45j . In 

total, the single-shot work associated to Steps b and c of 
the process is = ksT In2 {D^-^{pi\\t^) — 

with failure probability at most 2e —« 
2e, when e is small. 

To show consistency we now consider the average ex¬ 
pected work extracted per copy if the single-shot protocol 
is carried out on V —>■ oo i.i.d. copies of the system. In 
such a calculation the work computed is an average value 
which is why the average work contribution of the 

basis rotation in Step a, can be taken into account too. 
One obtains a total average work per copy of 

(IF) = (If(“)) -f keT In 2 lim lim — 

e—J-0 N—>-oo N 

= (IF(“)) + keT ln2 [i^(pi||r") - Dip^Wr^)] 

= (IF(“))+tr[pi H]+\nZ-kBTS{pi) 

- ti[p^ H]-\nZ + kgT S{p^) 

= kBTiSip^)-S{p)) 

= (W^opt), 

(F3) 

where we have used the quantum asymptotic equiparti- 
tion theorem for relative entropies gSillT] in the second 
line. i^(j|-) is the standard quantum relative entropy 
defined by D{p^\\t^) = tr[7y^(log77'^ — logr^)] and 
likewise for pi, where log is the logarithm to base 2. The 
quantities and their regularized ver¬ 

sion, the standard quantum relative entropy D, can be 
seen as different measures characterizing the distance be¬ 
tween two states. When applied here, they measure the 
‘distance’ between the thermal state and another di¬ 
agonal state in such a way that the operational meaning 
of this distance is given by the work one has to invest or 
is able to extract when transforming one into the other. 

The derivation shows that in the asymptotic limit the 
optimal average work is recovered from the single-shot 
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components. But it is important to realise that from 
Eq. (F3) one cannot conclude that the above single-shot 
process forming from pi is optimal. Going via the 
thermal state is just one option which is particularly con¬ 
venient in this case as the processes of maximal work 
extraction and work of formation from the thermal state 
have been treated in the single-shot scenario [SD] . It is an 
open question whether there are better single-shot pro¬ 
tocols for general thermodynamic transformations, see 
Fig.i^ & 1^. The introduction of “catalysts” in single¬ 
shot thermodynamics [5] provides a promising avenue to 
establish bounds on the work that can be drawn from a 
state with coherences during a projection in the single¬ 
shot setting. 


After making public our results on average work asso¬ 
ciated with removing coherences in thermodynamic pro¬ 
jection processes very recently a paper appeared [31] that 
derives the work that can be extracted when removing co¬ 
herences in a single-shot setting. In this paper the previ¬ 
ously mentioned framework describing the catalytic role 
of coherence in thermodynamics by Aberg |12j is used 
together with insights from reference frames in quantum 
information theory. These results are in agreement with 
our findings and strengthen our conclusion that coher¬ 
ences are a fundamental feature distinguishing quantum 
from classical thermodynamics. 


Appendix G: Quantum work fluctuation relation 


The average exponentiated work then becomes 




where „ are the transition probabilities for energy 

jumps starting in and ending in at time r. 

These probabilities are given by 


Pl^,n = tr[nL") En(0) po ni") yt j 


(G3) 


= y^'j, 

simplifying the exponentiated average work to 

(e^^) = Y yy^]. 

m n 

(G4) 

The completeness of the projectors, = 1, now 

finally results in the well-known quantum Jarzynski work 
relation 






(G5) 


where AF = is the difference of the equilib¬ 

rium free energies corresponding to the final and initial 

Hamiltonians, i.e. = —ksT In Sim¬ 

ilarly, the average work extracted from the system is the 
average energy difference between po and pr 


A common route of deriving the quantum Jarzynski 
equation is as follows 011112]. A quantum system is ini¬ 
tialised in a thermal state po = ’) 

for a given Hamiltonian with 

energy eigenvector projectors = |en°^)(ei°^|, at 

given inverse temperature j3 = l/(fcBr). Here F^p'^ = 
—ksT In is the initial free energy associ¬ 

ated with the initial Hamiltonian . The aim is to cal¬ 
culate the average exponentiated work, that the 

quantum system will exchange when undergoing a uni¬ 
tary process V that is generated by varying the Hamil¬ 
tonian in time, i.e. F[^*\ from to a final = 

The final state after the unitary is the 
non-equilibrium state Pr \= V pqV\ see Fig. ir. 

To identify the work for an individual run of the ex¬ 
periment the energy of the system is measured at the 
beginning, by projecting into and at the end, by 

projecting in the final energy basis Hm^ = |em ^)( Cm I. 
The (extracted) fluctuating work W identified with each 
transition —>■ |em^) is the (negative) observed fluc¬ 

tuating energy difference of the system 

AF = Fi”) - f(°) = -ly. (Gl) 


m,n 

= -^FMp);, + ^F(0)p0 

m n 

= -tr[FM p.]+tr[F(°) po] 

= —UiPr) + U{po) 


where p^^ := En Pm,n = ^APt Hm^ ] are the probabili¬ 
ties to find energies Fm^ when a measurement of 
is performed on pr- By construction, the initial proba¬ 
bilities to find energies En'^ are just the thermal prob¬ 


abilities, := EmPEn 
(^-PiE^ff -eE ) , 




The derivations of the average work, = 

U{po) — U{pr), as well as the average exponentiated work, 
i.e. the quantum Jarzynski equality, (e^^) = e~^^^^ are 
based on Eq. ( |G1[ ) which was made assuming that the 
initial and final state of the process that is being char¬ 
acterised are po and pr- There is no question that the 
mathematical details of the derivations of the above rela¬ 
tion are sound. Experimentally, there is however a need 
to acquire knowledge of the fluctuating energy to quan¬ 
tify the work and this requires the implementation of the 
second measurement, see Fig. Only after the mea¬ 
surement has been made can theoretical predictions be 
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tested. The measurement is an unavoidable non-unitary 
component of the overall experimental process. Specifi¬ 
cally, the ensemble state after the unitary, p^, is further 
altered by the measurement to result in the final state p-r, 
i.e. it is the state pr with any coherences in the energy 
basis of removed. 

While the experimentally observed average energy dif¬ 
ference is not affected by the measurement step, i.e. 
U{f]r) — U{po) = U{pr) — U{po), the entropy difference 
does change, i.e. S{f]r) — S{po) ^ S{pr) — S{pq) = 0. 
This means that the system may absorb heat, (Q®'*’®), 
during the measurement step, indicated in Fig. Wp- Its ac¬ 
tual value depends on how the measurement is conducted 
with the optimal heat positive, (QopD = kBT{S{r]r) — 
S{Pt)) > 0. Since AU = — (W) (Tl) this implies 

that in an experimental implementation of the Jarzyn- 
ski relation the work done by the system on average can 
be more than previously thought, with the optimal value 
being (VFopt) = + ksT - S{p,)). In 

the special case that the average heat is z ero it 

is possible (although not necessary) that Eq. (GI), and 
thus the standard Jarzynski expression (e^^) = , 

are correct. In particular this applies to classical mea¬ 
surements. We conclude that the suitability of identi¬ 
fying W = —AE, and hence the validity of the quan¬ 
tum Jarzynski work relation depends on the details of 
the physical process that implements the measurement. 

Quantum work fluctuation relations that have only one 
measurement [351135], instead of the two discussed above, 
offer a feasible route of measuring work fluctuations ex¬ 
perimentally. Instead of measuring separately the initial 
and hnal fluctuating energies, En'^ and Em\ to estab¬ 
lish their joint probabilities, this method acquires only 
knowledge of the joint probabilities by measuring energy 
differences AE directly. But also here is one final mea¬ 
surement, in general on a non-diagonal state, needed. 


Appendix H: Access to correlated auxiliary systems 

Similarly to erasure with a correlated memory m one 
can consider projections on a system S that is correlated 
with an ancilla A the experimenter has access to. Assum¬ 
ing a total Hamiltonian E[^^ = 0 1"^ -|- 1'^ 0 we 

denote the global initial state by p^^ and its marginals 
on S and A hy p^ = and p^ = tr^lp'^"^], 

respectively. 

A note on notation. For clarity we employ a 
slightly different notation here. The roles of initial state 
p and final state rj are the same as in the main text and 
the previous sections of the Appendix. However, now 
the superscripts of the final state p no longer denote the 
projection basis but the system for which p describes 
the state. For instance, p^ denotes the reduced state 
after the projection on system S alone. The same holds 
for the superscript of the initial state, p^"^,p^ and 
p^, and the Hamiltonians and H^. Only 


the superscript V of the mutually orthogonal rank-1 
projectors {H^j^ acting on system S is kept to indicate 
which basis is being projected in. 

For an initial global state p^^ of system and ancilla a 
local projection map on S results in a new global state 

p^^ = 1^) p^^ (H^01^). (HI) 

k 

Due to the properties of the projectors the marginal state 
on A is unchanged, 

p^ = trs[77^^] = ^ trs[(n^ (g) 1^) (H^ (g 1^)] 

k 

= ^trs[/^ (H^ (g 1^)] = trs[/^] = /. 

k 

(H2) 


The reduced state of the system becomes p^ = 
trA[p^^] = J2k^kP^^k = J2kPkPlk where pk = 
tr[/9'^n](’], and the conditional states on A after the pro¬ 
cess are denoted p^ = Pk^ trs[ (H^ (g 1^) p^^ (H^ (g 1"^) ] 
for all k. The global entropy change associated with the 
local projection is 

= S{p^^) - S{p^^) 

= S{{p,})+J2pkS{pt)-S{p^^) 

k 

= S{p^) - 5 (/) + S{p^) + J2pkS{p^) - S{p^^) 

k 

= AS'^ +6'^{A : S). 

(H3) 

In the second equality it was used that p^"^ is a classical- 
quantum-state and S{{pk}) = —J2kPk^kipk stands for 
the classical Shannon entropy [50] which is equal to the 
von Neumann entropy of p^ because the final state on S 
is a classical mixture of states from the projective basis. 
Here we defined a measure of correlations between the 
ancilla and the system, 5^{A : S) = S{p^) — S{p^"^) + 
'^kPk^iPk)! related to the quantum discord. It depends 
on the projectors {H^}*, and is always positive [33 ED- 
Thus the entropy change of SA can be bigger than the lo¬ 
cal entropy change, AS^ = S{p^) — S{p^), on the system 
alone. 

As shown before, Eq. ([^, the optimal extractable work 
in a thermodynamic projection process on system S alone 
is (Wopt) = ksT AS^ — AU^, where AS^ is the entropy 
change of the system and AU^ its change in internal en¬ 
ergy. This result stays intact when generalizing to pro¬ 
jections in the presence of ancillary systems if one takes 
the total changes of these quantities on 5'A instead of the 
change on S only. In the global process the total inter¬ 
nal energy change is equal to the energy change of the 
system only as the local state of the ancilla is unchanged 
and the total Hamiltonian is the sum of local Hamilto¬ 
nians. Thus using side information the overall optimal 
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extractable work amounts to 

(Wopt) = fcsT - AC/^ 

= fcsT AS'^ + keT S'^ {A : S) - AU'^ (H4) 
= {W,pt)+kBT6'PiA:S), 

where (kkopt) is the work of an optimal thermodynamic 
projection process without access to correlated systems, 

Eq. (||). 

Discord was first discussed in a thermodynamic con¬ 
text by Zurek m, where he related it to the advantage 
a quantum Maxwell demon could have over a classical 
one. In general the quantum discord, S{A : S), is defined 
as the minimum of S^{A : S) over all sets of projectors 
{n^}fc whereas in our case this set is fixed (see e.g. Modi 
et al. [S2] for a review). Therefore it is found that even 
for states with no quantum discord, usually referred to 
as classically correlated states, a difference in work asso¬ 
ciated with thermodynamic projection processes can be 
observed. This contrasts with the erasure process [T3] 
where an advantage could only be gained for highly en¬ 
tangled states. 

One may ask what global states on SA maximize 
(Wopt) for a given state on S. Expectedly, it can 
be shown that purifications of p'® yield the best improve¬ 
ment in terms of extracted work. Given p^ = ai\l){l\ 
any purification is, up to isometries on the purifying sys¬ 
tem HU, equivalent to = X); for some 

orthonormal basis of A. For such a state the con¬ 

ditional states on A after the projection, p^, are pure 
for all k which implies that they have zero entropy. This 
implies 

S-^iA : 5) = Sip^) - ^(|4-)(vl/|) +Y,PkSivi) 

(H5) 

= 5(p^)-0 + ^p,.-0 = 5(p^). 

k 

The optimal total extracted work from a purified state on 
SA in a thermodynamic projection process is therefore 
(Wopt) = ksT Sirj^) — AU'^ which can be shown to be 


the maximum for fixed p^ and projectors {11^}^. One 
way to see this is the following Lemma [53]. 

Lemma 1. Let p^^ be an arbitrary state on a bipartite 
system SA, and let {n^}fc be a complete set of rank-1 or¬ 
thogonal projectors on S. Consider the global state after 
the projection, rj^^ = X(fc(nfe' ® 1'^) p'^"^ (n? (g) l'^) = 
® where pk = tr[(n^ ® l^)p®^] are the 
probabilities to measure k on S and = Pk^ ® 

lA^pSA (n^(g)l^)] 

are the conditional states on A. Then 
S{p^^) > Y^PkSirit). (H6) 

k 

Proof. We model the process on S as an isometry 
\ijjk)^ <8) nf, where 5 is a copy of S and 
{\tpk)^}k is an orthonormal basis of S. The state af¬ 
ter applying the isometry is denoted 
and we note that tr^jp'^'^^] = Furthermore, isome¬ 
tries do not change (von Neumann) entropy and thus, 
S{p^^"^) = S{p^^). In addition, by construction of 
S the marginals of the final state on S and S have 
the same entropy: S{p^) = S{p^). Since S{p^^^) > 
|5'(,yS'^) _ S'(p'5)| > S{p^^) — S{p^) (see e.g. Nielsen & 
Chuang [T7|). Thus 

Sip^^) = S{r]^^^) > Sip^^) - Sip^) 

= Sip^^)-Sip^)=J2PkS{vt), 

k 

where in the last equality we made use of the fact that 
p^^ is a classical-quantum state. □ 

Going back to Eq. ( |H3| ) and applying the the above 
Lemma we see that in general (Wopt) = ksT AS^ — 
AU^ = kBT[S{p^) + j:kPkS{vi)-S{p^^)] - 
AU'^ < kBT S{p^) — AU'^, which proves that purifi¬ 
cations on SA yield the maximally possible extracted 
work. 
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